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ABSTRACT 

A complete description of dynamics of compact locally homogeneous universes is given, 
which, in particular, includes explicit calculations of TeichmiiUer deformations and careful 
counting of djTiamical degrees of freedom. We regard each of the universes as a simply 
connected four dimensional spacetime with identifications by the action of a discrete sub- 
group of the isometry group. We then reduce the identifications defined by the spacetime 
isometrics to ones in a homogeneous section, and find a condition that such spatial identi- 
fications must satisfy. This is essential for explicit construction of compact homogenoeus 
universes. Some examples are demonstrated for Bianchi II, VIq, VIIq, and I universal 
covers. 



1 Introduction 



In relativistic and observational cosmology, we often use a simplified spacetime model 
having restricted dynamical degrees of freedom. In particular, the well-known homoge- 
neous and isotropic (FRW) models in which the spatial sections are assumed to be 
homogeneous and isotropic, have been successful. On the other hand, a wider class of 
models, known as the Bianchi homogeneous models P, |^, Q , in which the spatial sections 
are assumed to be homogeneous but not isotropic, have been largely used in relativity and 
quantum cosmology 0, |5|, ^. In the models except type IX, each spatial section has been 
regarded as open. The open topology, however, is not a sole possibility. For example, 
the "open" model in the homogeneous and isotropic models, which has constant negative 
spatial curvature, and belongs to Bianchi type V, can be regarded to be spatially compact 
if spatial points are appropriately identified with some discrete subgroup of the isome- 
try group. The purpose of this paper is, in fact, to investigate a class of the models in 
which each spatial section has compact topology. The compactness of space is physically 
reasonable due to its finite spatial volume. 

The crucial point of the arguments in this paper is that the compactness of locally 
homogeneous space, in general, brings about new degrees of freedom of deformations, 
known in mathematics as Teichmiiller deformations. They preserve the local geometry 
but change the global one. This can be easily understood if we regard the Teichmiiller 
deformations of a compact locally homogeneous space as a homogeneous space (i.e., a cov- 
ering space) with varying identifications. A space spanned by independent Teichmiiller 
deformations, and its coordinates are referred to as the Teichmiiller space and the Te- 
ichmiiller parameters, respectively. We shall shed light on these degrees of freedom of the 
Teichmiiller deformations, which in fact have been often disregarded so far. The field of 
(2+l)-gravity was exceptional, but we shall take a somewhat different approach. The 
Teichmiiller deformations are of great interest, because they would carry part of the dy- 
namical degrees of freedom. We shall give a complete framework to do a concrete analysis 
concerning the Teichmiiller deformations. We also carefully count the total dynamical 
degrees of freedom. 
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In Ref.pi, referred to as I hereafter, we presented a treatment of three-dimensional 
compact homogeneous Riemannian manifolds, which will be a basis of our arguments in 
the present work. We, there, (1) gave the possible eight types (a~h) of homogeneous 
universal covers, which are closely related to Thurston's eight geometries[|, |10], (2) clas- 
sified compact quotients (al/1, b/1, etc.), and (3) gave Teichmiiller spaces by explicitly 
finding embeddings of covering groups in the isometry groups of the universal covers, 
which enable us to perform explicit calculations. To investigate the dynamics of compact 
homogeneous universes, we must first show how we can adapt such knowledge of compact 
homogeneous three-manifolds in the context of relativity in four dimensions. 

Our strategy to this will be as follows. We begin with considering a four- dimensional 
universal cover (^^^M,gab), which is a simply connected Lorentzian manifold, and then 
take identifications in {^^^M, gab) so as to make each 3-surface (Mt, hab) compact. To 
utilize our knowledge about compact homogeneous 3-manifolds, we translate the identi- 
fications in gafe) into those in {Mt,hab)- We find that the identifications must be 
extendible isometries of {Mt,hab), which have natural extention in {^^^M,gab)- We make 
(Mj, hab) compact by the action of a discrete subgroup F of the group EsomM^ C IsomM^ 
of extendible isometries, where IsomM^ is the isometry group of {Mt,hab)- Once given 
r C EsomMj on (Mt,hab), all necessary identifications in (^^^M,gab) are automatically 
determined by the natural extension of T. The quotient C^^^M, gab) is a solution of a lo- 
cal equation, e.g., Einstein's equation, if and only if {^'^^M,gab) is a solution of the same 
equation. 

The organization of this paper is as follows. In section H, we first give the defini- 
tion of Teichmiiller deformations, and then briefiy review the classification of compact 
homogeneous 3-manifolds given in I. In section ^, we establish the prescription for identi- 
fications, and discuss possible four-dimensional universal covers. We also discuss how we 
can eliminate the "gauge" degrees of freedom, and thereby we give how to find the dy- 
namical degrees of freedom. In section ^ we apply the framework of the previous section 
to concrete models. There, we give the time-development of the Teichmiiller parameters, 
establish the dynamical variables and give the number of the dynamical degrees of free- 
dom for each case. The final section is devoted to conclusions. We employ the abstract 
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index notation (See e.g. [Q) throughout this paper. 



2 Preliminaries 

We give definitions of Teichmiiller deformations in the first subsection. In the next sub- 
section, we briefiy sketch the classification scheme, given in I, of compact homogeneous 
3-manifolds, though we will not duplicate the results of the classification. This will be 
helpful for the subsequent discussions. For the explicit results, see tables 1 and 2 [|TT] in 



I, and section V of I. We consider only complete Riemannian manifolds, and shall drop 
the word "complete" hereafter. 

2.1 Teichmiiller deformations 

We define Teichmiiller deformations of a Riemannian manifold (M, hat) as follows. 

Definition 1 (Teichmiiller deformations) Let (M, hab) be a Riemannian manifold. 
Then, smooth and non-isometric deformations of hab ff^e called Teichmiiller deformations 
if they leave the universal cover (M , hab) globally conformally isometric. 

In this definition, a globally conformal isometry means a conformal isometry with constant 



conformal factor. For definition of coverings, see, e.g., Ref. ||12 . 

In other words, Teichmiiller deformations are deformations induced by variations of a 
covering group F of the universal cover (M, hab)- Here, a covering group F is a represen- 
tation (or an embedding) of the fundamental group ttiM into IsomM, the isometry group 
of (M, hab), where M can be realised as M/F. We denote the space of all covering groups 
as Rep (vTiM, IsomM). Note that not all variations of F correspond to independent Te- 
ichmiiller deformations. In fact, two Riemannian manifolds M/F' and M/F are isometric 
if 

r = aoToa-^ (1) 

holds for an isometry a G IsomM. F' is called the conjugation of F by a. If we denote 
the equivalence relation by conjugation as ~, the Teichmiiller space, Teich.{M , hab) , for 
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M and hat, is defined as 

Teich(M, hab) = Rep(7riM, IsomM) / ~ . (2) 

Although in this paper we are interested only in locally homogeneous metrics, it is 
worth noting that the definition of Teichmiiller deformations here concerns wider classes of 
metrics. For example, even though the universal cover admits only one Killing vector, the 
quotients can smoothly deform if we smoothly vary the identifications along the Killing 
orbit. 

2.2 The classification of compact homogeneous 3-manifolds 

We briefiy sketch the classification scheme, given in I, of compact homogeneous 3-manifolds. 

Consider a pair (M, G) of a manifold M and a group G acting transitively on M 
with compact isotropy subgroup. Note that we can construct a homogeneous manifold 
(M, hab) by first giving an arbitrary metric at a point p in M, averaging it by the isotropy 
subgroup, then finally sending it by the actions of G. Since G is transitive, the resulting 
metric hab is guaranteed to be homogeneous. The isometry group of (M, hab) would be 
isomorphic to G, or contain G as a subgroup of the isometry group. Note that if we 
give various metrics at p, then we obtain many homogeneous metrics on M. Hence the 
pair (M, G) can be considered as an equivalence class of homogeneous manifolds whose 
isometry groups are isomorphic to G, or contain G as a subgroup of the isometry groups. 
Such a pair (M, G) is called a geometry. If two geometries (M, G) and (M, G') have an 
inclusion relation G C G', then {M,G) is called a subgeometry of {M,G'). If geometry 
(M, G) is not a subgeometry of any geometry, then (M, G) is called a maximal geometry, 
and if geometry (M, G) does not have any subgeometry, then (M, G) is called a minimal 
geometry. 

Our starting point of classification of compact homogeneous manifolds is following 
Thurston's theorem [Q. 

Theorem 1 Any maximal, simply connected three-dimensional geometry which admits a 
compact quotient is equivalent to the geometry (M, IsomM) where M is one of , H^, 
S^, xR, H^x R, SL{2,R), Nil, and Sol. 
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A brief proof and accounts of the eight geometries in the theorem are found in Ref. |T3 • 

Note that this theorem concerns only maximal geometries. This seems, however, too 
restricted for our purpose. For example, while we are interested in the closed FRW 
universe, of which a spatial section corresponds to maximal geometry (S'^,S0(4)), we 
are also, and maybe more, interested in the Bianchi IX universe, of which a spatial 
section is subgeometry (5*^,811(2)) of geometry (S'^,S0(4)). Hence we should concern 
all non-maximal geometries which admit compact quotients, too. The above theorem, 
however, is still of great use, because of the following fact. That is, the "maximized" 
geometry of any non-maximal geometry admitting a compact quotient must admit a 
compact quotient, because the group of the maximal geometry must contain the covering 
group which makes the non-maximal geometry compact, and it must make the maximal 
geometry compact. This simple fact and Thurston's theorem lead us to the investigation 
of all possible subgeometries of Thurston's eight maximal geometries. Of course, we must 
note that not all subgeometries of the eight geometries admit compact quotients. We 
need to check explicitly that the subgeometry certainly admits a compact quotient. 

Although, as we have seen, the use of "geometry" is convenient to carry out classifi- 
cation of homogeneous manifolds, it is useful to switch to the conventional "Riemannian 
manifold" representation for physical applications. As we have noted, a geometry is an 
equivalence class of homogeneous Riemannian manifolds. Let R be the set of all homo- 
geneous manifolds belonging to a geometry, and let R be the quotient set of R by all 
possible diffeomorphisms and globally conformal transformations. Our classification for 
all possible universal covers admitting compact quotients is described in terms of ^'s. 
In I, we labeled those ^'s as "al, a2, b, c, h". For each element of quotient R, we 
choose a representative element (M, hah), and call metric hab the representative metric or 
the standard metric, whose explicit form for each type is also given in I. 

In getting the standard metrics, the Bianchi-Kantowski-Sachs-Nariai (BKSN) classifi- 
cation 0, 11^, |T3|] is useful, which concerns all minimal 3-geometries; 



Theorem 2 Any minimal, simply connected three-dimensional geometry is equivalent 
to {M,G) where M = R^, G =one of Bianchi I to Bianchi VIII groups; M = , 
G = Bianchi IX group; or M = x R, G = 50(3) x R. 



This is useful in that the invariant metrics for these groups are well known. These ge- 
ometries are called the BKSN minimal geometries. We take this opportunity to show 
the correspondence of Thurston's eight geometries and the BKSN minimal geometries, 



which is shown in table 1. (Such a correspondence was first pointed out by Fagundes 
without referring to subgeometry. ) 



Thurston's maximal geometries 


BKSN minimal geometries 


(E^JsomE^) 


(R3,B1), (R3,BVllo) 


{H^, IsomH^) 


(R3,BV), (R^BVII,) 


{S^, IsomS^) 


(53, BIX) 


{S^ X R, IsomS^ X R) 


X R, KSN) 


{H^ X R, IsomH^ x R) 


(R3,B111) 


(SL(2,R),IsomSL(2,R)) 


(R3,BVIII), (R3,B111) 


{Nil,lsomNil) 


(R3, BlI) 


{Sol, IsomSol) 


(R3,BVlo) 



Table 1: The correspondence between Thurston's maximal geometries and BKSN minimal 
geometries. BI is an abbreviation for Bianchi I group, and similar for BlI, Bill, etc. KSN 
is equivalent to IsomS^ x IsomR. The missing Bianchi types IV and Vl^ do not admit 
compact quotients. 

It should be noted that an R does not always correspond to a single geometry, though 
the converse is true. For example, consider Bianchi 1 minimal geometry (R^,R^), where 
the left R^ stands for the manifold homeomorphic to R^, while the right R^ stands for the 
3-dimensional translation group acting on the manifold R^. Then the corresponding R is 
the set of Riemannian manifolds (R^, /iafe)'s, where hat = hf^^{dx^)a{dx'^)h with h^^, being 
positive definite symmetric 3x3 matrix. All such {H^,hab)^s, however, are isometric to 
(R^,?7afe), where rjab is the standard Euchd metric rjab = 6^^{dx^)a{dx''')b. We thus have 
R = {(R^ , r]ab)} ■ On the other hand, it is manifest that maximal geometry (R^,10(3)) 
also gives rise to the same R, since 10(3) is the isometrics of rjab- (The R of this example 
is classified to type a2 in our classification. Similarly, Bianchi types II, VIIq, and VIq give 



rise, respectively, to universal cover types b, al, and f, on which compact models will be 
discussed specifically in subsequent sections.) 

Let us summarize the actual procedure for our classification. First, list up Thurston's 
eight geometries and all their subgeometries. Enumerate all possible compact quotients 
of Thurston's eight geometries, and check whether the subgeometries admit them. Find 
out ^'s for the subgeometries which admit a compact quotient. Note that this includes 
the explicit determinations of the standard metrics. This completes the classification of 
the universal covers which admit a compact quotient. The classification of the compact 
quotients that those universal covers admit can be carried out by using the results in 
Ref. [p!7| |T^. The Teichmiiller parameters are usually defined as components of the 
identification generators acting on the standard metric. Our classification, as a result, 
consists of the classification of universal covers, the classification of compact quotients, 
and explicit parametrizations of the Teichmiiller spaces. 



3 The framework of construction 

In the first subsection, we show a method of construction of compact homogeneous uni- 
verses, and then in the second subsection, we discuss the dynamical degrees of freedom of 
a system of compact homogeneous universes. Note, however, that the second subsection 
is not independent of the first one. Our construction of a system of compact homogeneous 
universes is completed by the discussions there. 

3.1 Universal covers and identifications 

We mean by a compact homogeneous universe a smooth Lorentzian 4-manifold {^^^M, gab) 
which admits a foliation by compact homogeneous spatial leaves (sections), and denote 
the universal cover of {^'^^M,gab) as {'^'^^M , gab) ■ It is important that {^'^^M,gab) inherits all 
the local properties from ('■^-'M, gafe)by a covering map. Hence, we can think that the local 
and the global geometries are carried by the universal cover C^^W, gab) and the covering 
map, respectively, and thus we can look into them separately. First, we shall consider how 
we can take the covering map, or "identifications", when a universal cover g^b) is 
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given. After that, we shall consider what universal covers are appropriate for our purpose. 

The identifications in gab) act on each homogeneous 3-section (Mj, hab) of t =const., 

where t parametrizes the homogeneous sections of (^^W, gab)- For simplicity, we omit the 
argument t of the metric hab as far as no confusions occur. Let IsomM^ be the isometry 
group of (Mj, hab)- It is very important to note that we cannot make, in general, the ho- 
mogeneous 3-manifold {Mt, hab) compact by the action of a discrete subgroup of IsomM^, 
since for (^^^M, ga^) to be a smooth Lorentzian manifold, the covering group, F, of the 
section needs to preserve the extrinsic curvature, as well as the spatial metric. To give a 
suitable prescription of compactification, we convert this requirement of the smoothness 
of ((^^M, gab) into the following statement; Since (^^^M, ga^) is obtained by taking iden- 
tifications in C^^W, gab), for any two points which are identified, there should exist an 
isometry of (^^^M, gab) {not of (M^, hab)) which maps one to the other. Hence, if we define 
the extendible isometry group of (Mj, hob)-, EsomM^ C IsomM^, as below, then we obtain 
a complete prescription for construction of a compact homogeneous universe, as shown 
subsequently; 

Definition 2 (Extendible isometry group) Let {Mt, hab) be a spatial section of {^'^^M , gab) ■ 
An extendible isometry is the restriction on Mt of an isometry of (^'^W, gab) which pre- 
serves Mf They form a subgroup of IsomMf. We call it the extendible isometry group, 
and denote it as Esom(Mt, ^^W), or simply EsomM^. Obviously, an extendible isometry 

a G EsomMt has the natural extension on '■'^-'M which is an element of Isom*^^™ and 
preserves Mt. We call such the natural extension on ^*^M the extended isometry of a, or 
simply the extension of a. 

Proposition 1 The identifications on an initial surface {Mt, hab) must be implemented 
in Esom(Mt,(4)M), 

F C Esom(Mt,(4)M), (3) 

to get a compact homogeneous universe out of a given four- dimensional universal cover 
{('^^M,gab)- Moreover, the identifications acting on whole {^*^M,gab) are determined by the 
action of the extension ofV on ^^^M . 
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For example, Kasner type metric, 

ds^ = -dt^ + t^P'dx^ + t^P'dy^ + t^P'dz^, (4) 

where pi ~ ps are constants, has Euchd spaces as its spatial part. Suppose (Mt, hat) is 
such a Euclid space, where the metric is given by 

df = fP'dx^ + e^'dy^ + e^-'dz^. (5) 

As long as a generic case is concerned where p^'s are all different, the (continuous) isome- 
trics of the 4-metric are only the translations generated by d/dx, d/dy, and d/dz, 
which form group R^. Since they all preserve the spatial sections, we find EsomoM^ ~ R^. 
(EsomoMf is the identity component of EsomMf , and similar for IsomoMf .) On the other 
hand, rotations generated by vector 

k, = -t-'^Pl-P^^y + tP^-P-X—, (6) 

OX oy 

and the similar vectors obtained by permutations of indices also form intrinsic isometry 
group, S0(3), of {Mt,hab)- This shows EsonioMt ^ IsomoMj ^ 1S0(3). On a homo- 
geneous section {Mt,hab), if we choose a covering group F in EsomMj such as to make 
(Mj, hab) compact, we obtain a compact homogeneous universe as a four-dimensional 
manifold through proposition |I|. 

We now discuss what universal covers are appropriate for our purpose. If the topology 
of ^^^M is fixed, then manifold is uniquely determined. Hence we only need to consider 
possible metrics, gat's. Let us consider the case where each homogeneous spatial section 
(Mf , hab) corresponds to a Bianchi minimal geometry. It is also straightforward to adapt 
the following description for the KSN minimal geometry. By the definition of our compact 
homogeneous universes, the metric should be of the form 

ds'^ = -N'^{t,x)dt^ + Kp{t){N'^{t,^)dt + a'^){N^{t,^)dt + a^), (7) 

where (t,x) are local coordinates, a" the invariant 1-forms, and a,(5, . . . run from 1 to 3. 
The spatial metric hai3{t)a°'a^ is, indeed, homogeneous on each section t =const.. 

For a metric with generic lapse and shift functions, the extendible isometry group 
EsomMj contains only the identity element so that we cannot use the prescription for 
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the compactification. One might argue that there would remain possibihty to find a 
discrete group T' for the compactification if the lapse and shift functions have some 
periodicity. However, such a discrete group does not contain continuous parameters. 
Since we are interested in the case of maximal number of continuous parameters in the 
initial identifications, we demand that EsomMf is transitive and therefore the lapse and 
shift functions are independent of the spatial coordinates x; 

ds' = -N\t)dt^ + V(t)(iV°(t)dt + (T°)(iV^(t)dt + af^). (8) 

This metric becomes the following form 

ds^ = -dt^ + h^p{t)a''a^, (9) 

by the induced map of a diffeomorphism which preserves each homogeneous section of 
t = const.. We shall mainly focus on this type of metrics hereafter. 

We here comment on an intuitive prescription for identifications, which is also useful 
particularly for the metric (^. Note that the normal geodesies emerging from a section 
Mt are uniquely defined, provided that they are parametrized by proper time r. We refer 
to the exponential map exp Tn°'{t) which is defined with respect to the normal vector field 
n''(t) on Mt as the normal map. (The image of Mt by the normal map exp™"(t) is not 
generally Mf for some t' when considering metric (|^. For metric (^, we of course have 
exp™"(t) : Mt —>■ Mt+T-) Obviously, if two points, a and b, on Mt are identified, any two 
points mapped by the normal map should continue to be identified, 

Vr G R; (exprn"(t))(a) ~ (exp rn"(t))(6). (10) 

Hence, we can determine how the identifications evolve in time, in terms of geodesies 
in a given four-dimensional universal cover. For the metric (y), since the hypersurface- 
orthogonal geodesies coincide with the t-axes, we immediately obtain the following useful 
proposition. 

Proposition 2 In terms of the coordinates [t, x) of metric if at the initial surface 
t = to an identification is specified as (to;x) ~ (to, ax), where a is a free action on the 
coordinate space, then at any time t we must have (t,x) ~ (t, ax). 
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That is, if we take identifications in EsoniMt on an initial surface of the metric (P), 
and describe them in terms of the spatial coordinates x, then the description of the 
identifications holds for any time t. 

By proposition one might conclude that no interesting global, i.e. Teichmiiller, 
deformations occur for the metric (^, since the identifications on each homogeneous 
section in terms of the spatial coordinates do not vary with time. However, this is not 
the case. Remember that the Teichmiiller deformations are defined with respect to the 
intrinsic geometries of the three dimensional sections (Mf, hat)- Variation of metric with 
time does cause Teichmiiller deformations with time, and if there exists difference between 
EsomMj and IsomM^, as in the example below proposition the situation becomes much 
richer. We will comment on this point again at the end of the next subsection. 

3.2 Dynamical degrees of freedom 

Let us consider a universe characterized by an initial data set {hab,Kab), where hab and 
Kab are the spatial metric and the extrinsic curvature of the initial spatial section M. 
To give a possible {hab, Kab) in the case that M has nontrivial topology and {hab, Kab) is 
locally homogeneous, we first need to cover M with some open patchs, define coordinates 
for each patch, and give a transformation function for each overlap of two patchs. That 
is, we need to set an atlas. After that, if we give an (/i^jy, K^j^^) on a point with respect to 
the coordinates defined in a patch, we can send {hf^u, K^^)^^^ to all points in the patch 
by the transitive group action, and finally assign the values of the data set to all patchs 
by the transformation functions and the group actions. Hence the information of the 
data {hab, Kab) is equivalent to the information of the value of (/i^^, K^^^) on a point and 
the way of taking transformation functions if patchs are fixed. However, it can be easily 
imagined that it is very difficult to count the number of possible {h^y^K^j^ys and the 
number of possible ways of taking transformation functions up to diffeomorphism. We 
can accomplish this counting, using coverings, as follows. 

As stated in the previous subsection, we think that a set, U , of universal covers, 
(*^^^M, gafe)'s, carries the degrees of freedom of local geometry like local curvatures, and 
the covering maps do the degrees of freedom of global geometry like Teichmiiller param- 
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eters. In this approach, it is evident that we need to ehminate the degrees of freedom 
of all possible diffeomorphisms in U. We introduce the equivalence relation in U by 
diffeomorphisms, and denote the resulting set of universal covers as U. If we fix a ho- 
mogeneous surface {Mt,hab) in a m = {^^^M,gab) e f/ ill, and suppose T C EsomMj 
makes (Mj, hat) compact, then we can identify a pair {u, T) with a compact homogeneous 
universe through proposition |I|. We denote the set of (u,r)'s of all possible F's for a 
fixed u as Cu- If F' = o F o holds for a G EsomMt, then the resulting compact 
homogeneous universe, (m,F'), is isometric to (m,F). In this sense, the freedom of taking 



conjugations of F by EsomM^ also corresponds to "gauge" freedom (cf. Sec. |2.1|) . Intro- 
ducing the equivalence relation in Cu by the conjugations, we get the quotient set Cu- 
Now, our dynamical system, C, of compact homogeneous universes is equivalent to the 
set C = {c\ c e Cu, u e tJ}. 

If we choose a homogeneous section arbitrarily for each element in C, we will have 
a set, /, of initial data sets, {hab, Kab)^s, on a fixed compact 3-manifold. In /, there are 
no elements which are isometric to each other, since for any different elements in C are 
non-isometric, and the development of an initial data set is unique. The set / is therefore 
equivalent to the set we considered at the beginning of the subsection. 

The following proposition is now trivial. 

Proposition 3 The number, dim (7, of degrees of freedom of a system of compact ho- 
mogeneous universes is the sum of the number, dim U , of the degrees of freedom of the 
four- dimensional universal covers up to isometry, and the number, dimC*u, of degrees of 
freedom of initial identifications, i.e. covering groups on an initial section, up to conju- 
gations taken by the extendible isometry group. 

Hereafter, by a construction of compact homogeneous universes, we mean a construction 
with explicit determination of representatives of the universes in the above sense, so that 
the number of arbitrary parameters in a universal cover should be dim U, and the number 
of arbitrary parameters in the identification generators should be dimCu- Note that for 
vacuum solutions of Einstein's equation for Bianchi class A [^(i.e. types I, II, VIq, VIIq, 
VIII, and IX) and type V, the metric components ha/sit) in Eq.(^ is "diagonalizable" , 
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i.e. becomes diagonal by diffeomorphisms. Hence we will begin with the diagonal form 
metric to construct compact homogeneous universes on Bianchi class A or type V. 

Note that we take conjugations for {Mt, hab) only by EsomM^ to obtain the initial 
identification parameters as stated in proposition ||, while the Teichmiiller parameters are 
defined with respect to conjugations by full IsomM^ (cf. SecpJ]). Roughly speaking, the 
difference between the freedom of EsomM^ and that of IsomM^ corresponds to the free- 
dom of giving initial 'velocities' of Teichmiiller parameters, as we will see more explicitly 
through the examples in the next section. 



4 Four compact homogeneous universe models 

In this section, we construct four explicit models of compact homogeneous universe, the 
b/1, fl/1, al/1, and a2/l models. For each model, we count the number of dynamical 
degrees of freedom and give the time-development of the Teichmiiller parameters. 

To get the Teichmiiller parameters of a compact section {Mt,hab), we need to com- 
pare two mathematical representations, i.e., {Mt, hab) with the covering group F, and the 
standard universal cover (M, h^^^) with the covering group. A, parametrized by the Te- 
ichmiiller parameters. (In I, the standard metrics were called the representative metrics.) 
F and A are generated by the same number, n, of generators, {gi} and {aj} {i = 1, ■ ■ ■ ,n), 
respectively, {g^} and {aj} satisfy the same multiplication rule of an extendible isometry 
group. We can get the Teichmiiller parameters by finding the automorphism of EsomM 
which relates the two sets of generators. We shall do this first for the b/1 model, where 
we will see the most typical calculation to get Teichmiiller parameters. Then the fl/1 and 
al/1 models follow. Finally, for the a2/l model, we present a different method in getting 
the time-development of the Teichmiiller parameters. 

Our universal cover metrics are synchronous (Eq.(PP) and diagonal (See Sec. [3.21). 
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4.1 The b/1 model: a compact model on Bianchi II geometry 

We start with the multiphcation rule of Nil (=Bianchi II group); 



/ 1 \ 

\3 I 



g' + h' 



\ 



J 



where g,h & Nil, and we shall use superscripts to denote the components of a group 

clement. We use the same components {x^,x^,x^) = {x,y,z) as coordinates of Mt. The 
action of Nil on is defined by the left action on {x, y, z) eNil. A Nil-invariant (diagonal) 
metric is given by 

dZ^ = /iiidx^ + /i22dy^ + /i33(d2 - xdy)^ (12) 



where {a 



3) are constants, i.e. independent of {x,y,z). The four dimensional 



universal cover metric of our concern is of the form 



(13) 



with haa being functions of t. The vacuum solution is, of course, known, but we proceed 
with calculation, leaving haa free, since they are complicated functions in the synchronous 
gauge and, moreover, it enable us to apply the result also to models other than the vacuum 
model. 

We consider manifold "b/1 {n = 1)", classified in I, which is probably the most 
stereotypical compact manifold modeled on Bianchi II geometry. The fundamental group 
TTi is given by (See Eq.(118) in I) 



= (^1, 52, 53; bi, g2]g^ \ bi, 53], [52, 53] 



(14) 



The topology of b/1 is illustrated in Fig.l. 

To represent the generators of tti, gfj's, in EsomM^ = Nil, we put them as 



g? 
\g? } 



3), 



(15) 
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Fig.l. The topology of b/1. Letters A and F show how the front and 
back sides are identified. The right and left sides, and the top and 
bottom sides are identified normally. The five arrows shows actions of 
QiS, illustrating the first relation in Eq.(f2), g^g2gi = 9192- 



and substitute these in the relations of vri (Eq.(p^). We then get the following; 



9i 



{91' J 



92 



92' 

y92' J 



93 



1 ^ 



^3" 



(16) 



where 5(3^ = gx^g2^ - g{^g2^ ^ 0. 



We then consider the possible conjugations by EsomMf = Nil. For the conjugation 
of (7j's hj h = {h^,h'^,h'^) G Nil (For typographical convenience, we sometimes write 
components of group horizontally), we have 



hgih ^ 



91 
91' 



9i' + h^9i' - h^9i^ 



, hg2h ^ 



92 

92' 



g2' + h'g2' - h'g2' 



hg^h ^ = 5-3. 



(17) 

We can make the third components of gi and g2 zero if we take h as = {gi^g2^ — 
91^92^)/ 93^1 = {gi^g2^ — gi^g2^)/g3^ ■ After all, our representation of tti in Nil reduces 
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to 



91 



V ° / 



92 



V " / 



93 



' ^ 



l9» / 



The nonvanishing four independent components in these (^j's determine the initial identi- 
fications in the universal cover with metric (|T3|) . 

To proceed further calculations, we here cite the definition given in I of the Teichmiiller 
parameters for b/1 and some related properties. We denote the standard universal cover 
as (R^, h^ah^i where the standard metric If^^ is given by (Eq.(75) in I) 



d/2 



dx^ + dy^ + (d^; — x&yf' . 



(19) 



Any compact homogeneous 3-manifold classified in b/1 is globally conformally isometric 
to manifold (R^, W^^^j where A is a covering group whose generators are given by 

/ n \ 

(20) 





\ ° / 



a2 



02 

V ° / 



as 








Then, the Teichmiiller parameters are r = (01^,02^,02^) (Eq.(129) in 1). We can see that 
the map 



S0 : 



y 



I 



Re 



\ 



(21) 



is a 1-parameter isometry for (R^, /i^ff ), where Re is the rotation matrix by angle 6, and 
(e is defined by 

Ce{x,y) = ^((x^ - y^) cos 9 - 2xysm9) sin 9. (22) 

We here remark that Se is not an element of EsomM^ but of IsomM^, and therefore 
EsomM^ 7^ IsomMf in the b/1 model. For an element h GNil, conjugation by se is given 
by 



Se 




h^ + Ce{h\h') 



(23) 
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Note that metric (p!2[) is rewritten as 

hnh22 



'33 



+ dy'^ + (d^' - x'dy'Y) 



(24) 



with 



hii f 
has 



(25) 



where h^a are regarded as constants. If we view this coordinate transformation as a 
diffeomorphism and drop the constant conformal factor of metric (123), the resulting metric 



coincides with the standard metric (|T^). This diffeomorphism is obviously an element of 
the HPDs p3| , H, from the form of metric (p4D, so that the transformation {x,y,z) —>■ 
{x' , y', z') is an (outer-) automorphism of Nil. The image of Qi^s, which acts on metric 
(H) (or metric (||)), is 



9i 




2/ 



( 

92 



92 



2/ 



92 

V ° / 



(26) 



Here, 



91 



1/ 



33 1 



-9i , 9i 



21 



22 



33 



-9l ,92 



1/ 



11 



' ^33 1 2/ 
7—92 , 92 
h22 



ihs 2 



(27) 



Generator is automatically determined by gi and g2 (see Eq.([T^)), so we will concentrate 
on gi, g2 and the images of them by automorphisms. Since Eq.(l26D is not of the form 
of Eq.(^), it does not yet give the Teichmiiller parameters. To get them, we take a 
conjugation of Eq.(P^ by the (full) isometry of Nil, which is given by Nil itself with sg. 
We can "rotate" the two-dimensional vectors {gi^',gi^') {i = 1,2) by conjugations by sg 
(Eq. (|23|)), leaving the third components zero by a conjugation by Nil like the way we 
obtained Eq.(|18|). So, we arrive at 

ai = hse^g[sg^h'^ = 



V 



{9i"y + {9i"y 
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and 



0-2 



where 



g2^' sin 6*1 + (72^' cos 61 




5'i 5'2 + 9i 92 



fi'l 92 



n 2'n 1' 



(29) 



cos 6*1 



fii 



1/ 



^{91"? + (fii^O 

and is an element of Nil. Using Eg. (p7| 
parameters; 



sin 6*1 



-9i 



21 



(30) 



V(fii'')' + (fii'')'' 
we obtain the final form of the Teichmiiller 



a2 



a2 



'^(fii^)^ + ^(fii^)^ 

'^22 'ill 



33 „ 1 „ 1 
— I 

22' 
/i33 



h 



~ 1 /^^l 32 + — fil fi2 



33 „ 2„ 2 
hll 



(31) 



Oi^ Vhnh22 

In Eq.(^), parameters 91^,91^,92^,92^ and hence (73^ are constants, and /laa's are 
functions of t. The metric components /iii,/i22, and /iss are determined by substituting 
Eqs.([TB|) and (|T2|) into Einstein's equation, and moreover we must exclude the degrees of 
freedom of HPDs from them (See Sec. 3^). Hence, the number of free parameters that 
the metric components can have coincides with the known number of degrees of freedom 
of the conventional (open) Bianchi models . For the vacuum Bianchi II, the number of 
free parameters in the metric functions is two, i.e., dimf/ = 2. With the four parameters 
specifying the initial identifications, dim Cu = 4, the total number of dynamical degrees 
of freedom of the present vacuum b/1 model is six. (cf. proposition |^) The dynamical 
variables are the Teichmiiller parameters ai^, a2^, and 02^, and the total volume 



V = (fi'3^)V^11^22/i33- 



(32) 



Remember that the Teichmiiller parameters are defined with respect to the standard 
universal cover which is isometric to the universal cover {Mt, hab) up to a global conformal 
factor. In fact, it is clear that, if we know the values of them, we can completely construct 
the original compact 3-manifold. 
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Additional remarks : 

We end this subsection with some common remarks to the subsequent subsections, 
which will be helpful for understanding the rather unique concept employed throughout 
this section. Since as in SecH any compact locally homogeneous 3-manifold can always be 
smoothly specified by some Teichmiiller parameters, curvature parameters, and the vol- 
ume, we are regarding them, rather than haaS, as the dynamical variables of the compact 
homogeneous universe. In the b/1 case, they are the four parameters (cii^, 02^, 02^, v), as 
pointed out. (There are no curvature parameters in this case). One remark we want to 
emphasize here is that we are on the standpoint that we do not ask whether or not such 
the dynamical variables directly fulfill some dynamical differential equations, since we do 
not need them to obtain the time-development of the "dynamical variables". (We will 
however discuss this problem in a separate work. See also Sec.^.) Of course, this is not to 
say we do not need Einstein's equation. The role of Einstein's equation in our calculation 
is to fix the universal cover, i.e. to fix haaS. 

One of the main conclusions in this subsection was that once the universal cover 
is set fixed by Einstein's equation, we automatically get the time-development of the 
dynamical variables through Eqs. ([3T|) and (p2D. In this sense, Eqs.(0) and ( |32D are the 
kinematical relation between the universal cover and the dynamical variables. Another 
remark is therefore the fact that even if some matter fields are included and the form 
of haaS accordingly vary, Eqs. (pT]) and ( p2D are invariant. Correspondingly, while the 
degrees of freedom, dim?7, of the universal cover may vary, those, dimCu, of the initial 
identifications are invariant. 



4.2 The fl/l(n) model: the compact Bianchi VIq model 

The multiplication rule of Sol (=Bianchi VIq group) is given by; 



9' 



\3 J 



7)3 



1 \ 



(33) 
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where g,h G Sol. We can easily check that the 1-forms 

1 , . o 1 



a 



e^dx + e ""dy), a = dz 



(34) 



V2' V2' 
are invariant under the left action of Sol if (x, y, z) is identified with an element of Sol. 
We therefore have the following invariant metric 



dl' = h,,{ay + h,^{a'Y + h,,{a'Y . 



(35) 



As in the case of Bianchi II, the four dimensional universal cover metric is Eq. ([13| ) with the 
above dP. In contrast to the Bianchi II case, the isometrics and the extendible isometrics 
coincide; IsomM^ = EsomM^ = (Sol plus three discrete elements). For future use, we 
present one of the three discrete elements here. It is 



h : {x,y,z) {-x, -y,z). 



(36) 



The fundamental group tti of a compact manifold modeled on Bianchi VIq geometry 
is given by (See Eq.(145) in I) 



TTl = (91,92, 93] [91,92], 939193^92^^ 939293 ^ 9192^' 



(37) 



where \n\ > 2. Because of the coincidence of the two isometry groups, we need not do 
new calculations to find embeddings of tti in EsomM^ other than those shown in I. We 
simply show the results. 
Let 



C3 = In 



n 



1 



(3^ 



and so are — e '^^ and 



-1 n 



If n > 2, then e '^^ and e^^ are the eigenvalues of matrix 

— e^^, if n < —2. Let {ui, vi) and {u2, V2) be the normalized eigenvectors corresponding to 
the two eigenvalues, i.e., 

/ 



Ml 
V ^1 



1 



V 



J 



IU2 

\ ^2 



1 



/ / n— \/??— 4 \ 
V -2 

V V — -2 — ) 



(39) 
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Then we can embed the generators of tti in EsomMt as (See Eqs.(156) and (157) in I) 



9i 



for > 2, or 



/ \ 

aQUi 

aoU2 




, 92 



aoU2 




aoV2 




, 93 



V 





Vcsy 



(40) 







(41) 



9i — aoU2 ? 92 — aoV2 ; 93 — h o 



for n < —2, where h is defined in Eq.(^). We thus find that the parameter for the initial 
identifications is only in Eq. (|40|) or Eq . (p|) . 

Before giving the time-development of the Teichmiiller parameters, we take this oppor- 
tunity to present a pictorial account to the topology of a compact Sol: f/l(?2). Manifold 
f/l(n) is a torus-bundle over S^. The relation (|37|) implies that gi and g2 generate the 
fiber torus. We can understand the topology of f/l(^) by observing the gluing map gen- 
erated by which maps generators of a torus to another generators of a torus. From the 
relation (BTI), we observe that 



9i = 939193 ^ 



92, 92 = 939293 ^ 



9291'. 



(42) 



This means that gs maps a parallelogram spanned by gi and g2 on an x-y plane to another 
parallelogram spanned by gi and g2 on another x-y plane (Fig. 2 (A) and (B)). If we "cut" 
and translate the second parallelogram by the actions of gi and g2, then we can take a 
fundamental region as a parallelopiped of which the bottom and top surfaces are spanned 
by gi and g2 (Fig. 2 (B) and (C)). When identifying the bottom surface to the top, the 
surface is stretched by n times, and then folded. (A geodesic congruence along 2;-axis in 
f/l(n) will therefore behave extremely chaotically after a journey over some periods.) 

Let us return to the operation to identify the Teichmiiller parameter. Note that we 
can transform the spatial metric (pSf) into 



h 



33 



(43) 
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Fig. 2: The topology of f/l(n = 4). gi and g2 span a parallelogram in a 
x-y ,say z = zq, plane (Fig. (A)). The stretched parallelogram, shown in 
Fig.(B), spanned by gi = 52 and ^2 = dtdi^ in z = zq + c^ plane can be 
identified by the actions of g\, 515^^, gi92^ with the parallelogram 
spanned by gi and 52 in a way that letters A~E in Fig.(B) label the 
corresponding regions. The corresponding regions in Fig. (A) are labeled 
by the same letters. So, the up and down sides of the parallclopipcd 
in Fig.(C), which shows a fundamental region of fl/l(4), are identified 
as indicated by circles and triangles. The front and back sides, and the 
right and left sides are identified in the trivial way. 
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where 



V 



(fellfe22)'/^ 

Z 



y 



(44) 



The transformation {x,y,z) {x',y',z') defined by Eq. (|4^) is an automorphism of SoL 
So, the appropriate action of identifications on {x',y',z') is given by 



9i 



( \ 



33 



aQU2 




, 92 



{hxih22Y" 



( \ 



33 



aoV2 




> 93 







(45) 



for n > 2, or (jfg = /io(0, 0, C3) for n < —2. This aheady coincides with the parametrization 
of the Teichmiiller space, i.e. the only Teichmiiller parameter is 



a{t) = j== — ao. 



33 



(46) 



It is worth noting that we can observe from Eq.(|45|) that the Teichmiiller deformations of 
fl/l(n) are the variations of the ratio of the area of the fiber torus to the length of the 
base 5*^ 

The dynamical variables (in configuration space) are the Teichmiiller parameter a, the 
3-volume v = /iii/i22^33(tto)^ up to a function depending on n, and the curvature control 
parameter A = ln(/iii//i22)- The number of dynamical degrees of freedom is four; one is 
for ao, and three is for the parameters contained in the four- dimensional universal cover. 



4.3 The al/1 model: a compact model on Bianchi VIIc 

The multiplication rule of Bianchi VIIq group is given by; 



/ .1 \ 



( 



V 




(47) 



where g,h & Bianchi VIIq group, and Rq3 is the rotation matrix by angle g^. The 1-forms 



= cos zdx + sin zdy, = — sin zdx + cos zdy, = dz 



(48) 
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are invariant under the left action of Bianchi VIIq. Hence the invariant metric is 



(49) 



As usual, the four dimensional universal cover metric is Eq.(|l3D with the above d/^. Since, 
as in the Bianchi VIq case in the previous subsection, the isometrics and the extendible 
isometrics of (Mj, hab) coincide, IsomMj = EsomM^ = (Bianchi VIIq group plus three 
discrete elements), our calculations to do will be similar to those in the Bianchi VIq case. 

The compact 3-manifold we consider here is al/1, which is homeomorphic to the 3- 
torus T^. The three generators, gi, g2 and g^, of the fundamental group of al/1 are hence 
all commutative. 

The embedding of the fundamental group in EsomM^ up to conjugacies by EsomMf (= 
IsomMj) is already given in I, which reads 



91 





2ln 



( 0^ ^ 

92 



92 



92^ 

2mTX 



^3 



9^^ 
2n'K 



(50) 



where l,m,n are integers. 

To obtain the time-development of the Teichmiiller parameters, we, as usual, first note 
that we can transform the spatial metric (H^) into 



d/2 = h 



where 



//ill 



Ih 



33 I \i j^{cosz'dx' + sinzdy')"^ + ^ j^{—smz'dx' + cosz'dy'Y + dz'^ ) , (51) 



V 



(fellfe22)'/^ 

Z 



(52) 



Here, the metric (|5lD coincides with the standard metric, given in I, of Bianchi VIIq up 
to global conformal factor. 

Since the transformation (x, y, z) (x', y\ z') is an automorphism of the Bianchi VIIq 
group, we can easily obtain the actions on the standard metric ( |5lD of the generators (|50D . 
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We immediately get 



ai = 



/ 1 \ 


y 2Z7r y 



, 02 = 



/ 1 \ 

2m7r 



, 03 = 



/ 1 \ 

as' 
2mr 



(53) 



where 



(h h "l-^/^ 

a/ = ^ = (1, 1), (2, 1), (2, 2), (3, 1), and (3, 2). 



(54) 
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Here, Oj's are the images of (^j's by the automorphism, and ai^ ~ a^^ are the Teichmiiller 
parameters. 

To summarize, the five constants ~ 53^ determine the initial identifications on an 
initial surface in (^^W, ga^)- The universal cover {^^^M,gab) have three arbitrary parame- 
ters in vacuum, and hence the number of the dynamical degrees of freedom is 8 (= 5 + 3) . 
We have seven dynamical variables; the five Teichmiiller parameters, the curvature control 
parameter A = ln(/iii//i22), and the 3-volume v = 1/122/^33 det(5fi, g2, gs). 



4.4 The a2/l model: a compact model on Bianchi I 

In this subsection, we give the time-development of the Teichmiiller parameters of the 
b2/l model (~ T^) on vacuum Bianchi I, by a method other than the one finding an 
automorphism explicitly. This is done by calculating invariants under the automorphisms, 
like lengths of minimal loops and angles between the loops of the compact homogeneous 

manifold. To this end, we introduce a matrix representing such invariants as follows. 

Let r be a covering group acting on {M,hab)- For a G F, p e M, let the map 
7a_j, : R ^ M be the geodesic satisfying 



7a,p(0) = P, 

7a,p(l) = a{p). 

Then, we define the map Vp : F — > by relating F to the geodesic generator at p; 

d7a,p(A) 



(55) 



Vp{a) = 



dA 



(56) 



A=0 
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Finally, let all the independent generators of F be Oj (i = 1, • • ■ , n). Then the loop matrix 
defined by 

Hij{p) = h{vp{ai),Vp{aj)) (57) 

will contain all the information concerning the global geometry of M = M/F. Here, 
we have dropped the abstract indices of the metric and the vectors in the r.h.s.. It is 
worth noting that the p-dependence of Hij{p) decides whether M is locally homogeneous 
or globally homogeneous. That is, if Hij{p) is independent of p, then M is globally 
homogeneous, and if not so, then homogeneity of M is local. 

We are in a position to calculate the time-development of the a2/l model in vacuum. 
Our four dimensional universal cover is the Kasner solution = (R^,gafe), where with 
the usual coordinates {t,x,y,z), gab is given by Eq.(|) with a = pi + p2 + Ps = 1 = 



(Pi)^ + {P2y + {PsY- Each homogeneous spatial section is given by (R^, /lat) with h. 



ab 



being Eq.(^). The covering group may be generated by three commuting generators, for 
which we write as 



9i 



J 



1 ~ 3). 



(58) 



Here, the action of gi on p = (xq, yo, Zq) on (R'^, hat) is given by 



9i 



^ 9i^ 



9i 



yo 



9i + Zq 



(59) 



Since the extendible isometry group of each slice is isomorphic to the commutative group 
R'^, we see that no nontrivial conjugation occurs. This implies that we cannot simplify 
the components of (/j's more than the original form of Eg. (|58|) . 

We can at this point count the dynamical degrees of freedom of the present model. 
Immediately can we see that the Kasner parameter carries the part of dynamical degrees 
of freedom in the universal cover, dimf/ = 1, while (yfj^'s in Eq.(|58|) carry the part in the 
initial identifications, dimCu = 9. The total dynamical degrees of freedom is therefore 
10 (= dimf/ + dim (7^). (cf. proposition |^) 
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Now, return back to the procedure to get the Teichmiiller parameters. We can easily 
find the geodesies satisfying Eg. (|55D , and then get the generator at p as follows. 



79„p(A) 



9i^^ + Zq 



9.' 



j 



(60) 



From this, we have 



= jlt'^'^9^''9r■ 



(61) 



On the other hand, any fiat {a2/l) can be implemented in the standard Euclid 
metric 



with six Teichmiiller parameters P3[ in three generators 



(62) 



ai 



\ " / 



^2 



V " / 



as 



as 



(63) 



Components of Vp{ai) is the same as Eq.(KD|) with replaced by a^-^ with = ai 



02 = 0. Using Eq. (|6^) as h, we have 



11 2 2 "^3 

^ (air 



(64) 



(a2^)^ + (a2^)^ 
y (sym.) {a3y + (03^)2 + (03^)2 J 

We set equal the two expressions Eqs. (|51D and ( pTD to get Oj-' as time functions with 
initial parameters t^i". Elementary calculations lead to the following results; 



lai 



(^2 



1\2 



2\2 



AV(aii)2, a,' = -a2\Y.t'^'''"'^93V)/^\ a,' = {detg)h'yA\ (65) 
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where = Ea ^^^''"^"^^3")^ (J = EaPa = 1, and is the (i, Q;)th cofactor of the 
matrix (g'i"). It would be useful, especially in getting 03^, to note that the determinant 
of Hij is given by 

det H = (det = {ai^f{a2^)\as^f. (66) 

5 Conclusions 

We have given a general method of construction of compact homogeneous universes. This 
is accomphshed by taking identifications in a universal cover {'^'^M,gab)- The universal 
cover must satisfy Einstein's equation, and the degrees of freedom of all the possible 
diffeomorphisms must be subtracted. The identifications in the universal cover are im- 
plemented by a discrete subgroup of the extendible isometrics, EsomMj. At this stage, 
one takes the conjugations by EsomM^, and finally we can obtain a system of compact 
homogeneous universes which is free from any diffeomorphisms, i.e., the free parame- 
ters in the metric and in the identifications are the dynamical degrees of freedom of the 
system. This method of construction is evidently applicable to any system of compact 
homogeneous universes (i.e., compact models on the Bianchi class A, class B, and the 
Kantowski-Sachs-Nariai models) . 

We have considered the dynamical variables of the system to be parameters specifying 
spatial sections completely. In this sense, the Teichmiiller parameters are dynamical 
variables, as well as the 3-volume and the possible curvature parameters. It is important 
to note that the number, /, of dynamical degrees of freedom is less than double the 
number, d, of dynamical variables. As we have seen in the explicit examples, not all of 
the initial velocities can be arbitrarily chosen. In some cases, / is less than 2d by 2. This 
could be explained by the Hamiltonian constraint. In the other cases, however, / is less 
than 2d — 2. These could be well understood if we study whether the dynamical system 
admits a canonical structure. This is also needed in canonical quantization of compact 
homogeneous universes. This problem will be discussed in a separate work. 

Although we focused on the time-developments of the Teichmiiller parameters and the 
dynamical degrees of freedom, our framework of compact homogeneous universes should 
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be useful in wide variety of problems in astrophysics, observational cosmology, fundamen- 
tal problems of relativity, quantum cosmology, and quantum gravity. For example, the 
behavior of geodesies in a compact universe can become chaotic (cf. Sec. [4.2|) , which fact 
may provide an interesting cosmological model. The problem of strong cosmic censorship 
for compact homogeneous universes is also of great interest, which is being investigated 
by some groups (See, e.g., PB|). 
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